Parametric Excitation of Alfven Waves by Gravitational Radiation 
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We consider the parametric excitation of Alfven waves by gravitational radiation propagating on a 
Minkowski background, parallel to an external magnetic field. As a starting point, standard ideal 
MHD equations incorporating the curvature of space-time has been derived. The growth rate of the 
Alfven waves has been calculated, using the normal mode approach. Various astrophysical applica- 
tions of our investigations are discussed, and finally we demonstrate that the coupling coefficients 
of the interacting modes fulfill the Manley-Rowe relations. 



I. INTRODUCTION 

Interaction of electromagnetic fields with gravitational 
radiation has been studied by several authors [1-19]. Be- 
sides a purely theoretical interest in such phenomena, 
there is a number of different applications. For example 
in astrophysics [1-5], in cosmology [6] and under labora- 
tory conditions [7-10], where - in the later case - the goal 
is to find suitable mechanisms to detect gravitational ra- 
diation. Furthermore, there are many examples of gravi- 
tational wave interaction that may take place in plasmas. 
This has been studied during the eighties by a group at 
the Kazan School of gravitation (see for example [1,10,11] 
and references therein) and more recently by Refs. 4, 12, 
13 and 14. 

In Ref. 12 it was shown that parametric excitation of 
high frequency plasma waves by gravitational radiation 
may take place. Due to the frequency matching condi- 
tions, however, the plasma must be very thin for that pro- 
cess to be possible, and the amount of energy transfer is 
therefore limited. In the present paper we will thus con- 
sider parametric excitation of low frequency MHD waves 
by gravitational waves, which - in contrast - may take 
place in a comparatively dense plasma. The relevance of 
this problem for the conversion of gravitational wave en- 
ergy to the plasma inside supernovas has previously been 
discussed by Ref. 14. However, due to the complexity of 
the physical situation, a highly idealized model will be 
studied, where a one dimensional monochromatic gravi- 
tational wave - superimposed on a flat background met- 
ric - propagates through a homogeneous two component 
plasma. 

The organization of the paper is as follows: In section 
II idealized MHD-equations incorporating the effects of 
the gravitational wave are derived, starting from covari- 
ant two-fluid equations. In section III parametric exci- 
tation of shear Alfven and magnetosonic waves are con- 
sidered, the three wave coupling coefficients are derived 
and the growth rate is found. By adding a phenomeno- 
logical resistivity to the equations, the threshold value of 



the gravitational amplitude is also calculated. Finally, 
in section IV, our results are summarized and theoreti- 
cal considerations like energy conservation properties and 
the fulfillment of Manley-Rowe relations as well as pos- 
sible applications are discussed. 



II. RELATIVISTIC MHD-EQUATIONS 

In order to obtain general relativistic fluid equations 
governing a plasma we begin by considering a system con- 
sisting of a charged perfect fluid and an electromagnetic 
field [20]. Introducing the restframe scalar quantities: 
mass density (or rather energy density times 1/c^) P(m)j 
charge density /9(q), pressure p and the 4- velocity field 
(or fluid velocity) m^ = dx^^/dr and 4-current density 



J' 



Pii)^ 



where t is the proper time and x^ coor- 



dinates in the lab frame, this system is characterized by 
having the energy- momentum tensor T^'^ = Th^^. +2^/lm) ' 
where 
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H) = (P(m) + ^) 
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rf"" , = — I Ff'^F^'' - -g^"'F''^F^ 
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and F^i, is the electromagnetic field tensor satisfying 



Maxwell's equations 



i^'^r. = -w^ 
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We have adopted the convention that Greek suffixes 
/i, v, .. have the range 0, 1, 2, 3 and i,j, .. have the range 
1,2,3 and the metric tensor g'"' has the signature (H 

-)■ 

The conservation laws of the system follows from that 
the 4-divergence of the energy-momentum tensor van- 
ishes, i.e. T^f = 0, and with the use of Maxwell's equa- 
tions one gets 



{P(ra)U^)■,^. + 4<M = 



(3) 



(m(i) + m(2))nu'^u'.^ = 5'^pj + F^^ 



(6) 



(P(r: 
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yunp.^ + F-'^j^ 
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where eq.(g) is obtained by projection along the 4- 
velocity u'^. This equation is identified as energy balance 
in the rest-frame of matter and gives the equation of con- 
tinuity (mass conservation) in the non relativistic limit. 
Equation (0) gives for j/ = the energy balance (modulo 
the content of (0), i.e., energy balance in the non rela- 
tivistic limit) and for v = 1,2,3 momentum balance. 

Under the conditions of low internal energy the fluid 
description of a plasma can be simplified by putting, for 
each particle species of the plasma, P(m) = mn where n is 
the restframe particle density and m is the particle mass. 
Also without these restrictions we may put yO(q) = qn 
where q is the particle charge. Suppose we have a plasma 
consisting of two species of particles oppositely charged 
(i.e. (71/(72 = ~l)j but, in general, with distinct masses. 
For each species we assign a fluid satisfying the equations 
(^ and (|j). The fluids are assumed interpenetrating and 
interacting through the electromagnetic field and, in gen- 
eral, the gravitational field. We neglect the effect of par- 
ticle collisions. If we assume non-relativistic pressure, i.e. 
such that mn ^ p/c^, and non-relativistic fluid veloci- 
ties - in the sense that we may neglect quadratic terms 
in 1/c in the z^ ^ components of eq.(0) - then we have, 
for each of the two fluids, equations for conservation of 
particles (or mass) and momentum in the form 



(nu^).^ = 



^P,3 + F\<inu'' 



Maxwell's equations remain the same if we let j^ be the 
total current density. 

Under the conditions that for both species dt ^ uj <^ uJc 
and C^ <C c^, this two-fluid description can be cast 
into a set of single-fluid equations. By dt ~ w we 
mean that a characteristic frequency, uj, can be assigned 
to the time variations in the dynamical quantities [21]. 
We use cOc = \q\B/m for the cyclotron frequency, B 
is the magnetic field intensity (which can be obtained 
from the Lorentz frame components of the electromag- 
netic field tensor) and Ca is the Alfven velocity, de- 
fined by C\ — B'^/^ouin. Under the above conditions 
it follows that the two fluid velocities are approximately 
equal and that the particle densities may be regarded 
as exactly equal. Furthermore, we let the equation of 
state be the one of isothermal compression and assume 
\u^u^.^\ ^ \F'^Qqnu°\/mn oc \F^,qnu^\/mn (meaning 
that the electric and magnetic forces approximately bal- 
ance each other). In light of this, we obtain the following 
set of single fluid equations (MHD-equations) 



(nu");^ = 



(5) 



F\u^ = 



nq(i) ■ 



m{2)P{l)^0 - W(i)P(2),j- 

"^(l) +"i(2) 



= kuTnA 



(7) 



(8) 



Equation (0) and (o) are obtained by adding the two 
particle conservation equations and the two equations of 
momentum balance, respectively, and setting the veloc- 
ities equal (when added), letting j be the total current 
density, p the total pressure and T = T^i) -l-T(2). The suf- 
fixes (1) and (2) refers to the two distinct particle species. 
By subtracting the two equations of momentum conser- 
vation, where the terms u^u*.„ have been neglected, one 
finds eq. (1^), which we refer to as the generalized Ohm's 
law. 

Note that if the right hand side of eq. (R) is negligi- 
ble as compared to F^f^vP and F^:U^ this equation simply 
reads F'^ u^^ = and we can then refer to the single- 
fluid equations as the ideal MHD-equations. In the fol- 
lowing wc limit ourselves to this case. Note that the 
MHD-equations are not independent of Maxwell's and 
Einstein's field equations, since we used Maxwell's equa- 
tions and T'^^ = 0, that follows from Einstein's equa- 
tions, in deriving them. 

We now consider gravitational radiation on a 
Minkowski background treating the plasma as a testfluid. 
Thus, the plasma back scattering effect on the gravita- 
tional fleld is lost. The gravitational radiation is chosen 
to be weak gravitational waves in the transverse traceless 
(TT) gauge propagating in the x'^-direction. This plane 
wave solution of the linearized Einstein field equations 
can be written 



ds' 



c^dt^ 



(1 - h+)dx^ 
2hxdxdy — dz 



(1 + h+)dy^ 



where h = /le'*'''^ -I- c.c. and \h\ ^ 1 with the wave 
vector [k^] — {uj/c, 0, 0, k) satisfying the dispersion re- 
lation k^kfj^ — 0. In all the following calculations we 
neglect terms that are quadratic in h or higher. The 
metric tensor can thus be written as g^^, = t]^^ + h^^ 
where 77^^ is the metric tensor of Minkowski space and 
hfj,i, represents the small, \h^^\ ^ 1, fluctuation in the 
gravitational field. 

The nonzero Christoffel symbols are then calculated to 
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where h = dh/d^ and £^ = x^ — x^. 

By expanding the covariant derivative, eq.(H ) becomes 
[nu^^),^ = as a result of F^^^ u"^ — and - noting that 

- eq.(g) reads 



G' = -'5>j-/i>j+0-'' 



(9) 



where we have introduced G* = mnT^j^u'^u'^ and tti = 
W(i) +m(2) . 

Next, we perform the same expansion in MaxweU's 
equations and rewrite them in terms of the electromag- 
netic field tensor in the form F'^. The idea is to express 
all field tensor terms in the same form, preferably the one 
that gives the most simple expressions. We can separate 
F'^f^ — ~^QJ^ into two equations. Setting /x = we 
obtain a Poisson-like equation which we discard - since 
in the MHD-regime j" « 0. Setting ^ — i we read off 
"Ampere's law", 

S^'F^^ = ^,of - {h^'F^)^, + r^^g^^Fl + T^^g^'^F^ 

(10) 

From eq.pp, which by symmetry in F'^'^ is equivalent to 
F^v,(j + -fl/cr,^ + Fcy^^u = 0, we obtain a number of trivial 
identities, a generalized equation for V-B and "Faraday's 
law": 



-^^0,2 ^ -^^2,3 + -^3,0 - (^x^3 - h+F 



-F\ 
F- 



0,1 



P 1,3 ~ -^3,0 - -{h+F^ + hxF^)^o 
{h+F\ + h^F%)^o 
+ {h^Fl~h+Fl),, 



^ 1,: 



F 



2,0 



For notational purposes it is convenient to introduce an 
abstract basis {x, y^ g}. The one-fluid equations and the 
Maxwell's equations above can then be written in a vec- 
tor representation with an algebraic structure identical 
to the Euclidean. We define 



X = a;x + yy + zz — X x + a; y + x z 
V = Uj^x + VyY + VzZ = ?i X + ti y + u z 



J = jx^ + iyY + izi- = j^x + py + fz 



A 



E = E,j:'k + Eyy + E,z = cF\x + cF\y + cF^^z 
B = B^St + By-y + B^z = F%x + F\y + F\z 



V = 9^x -I- 9^y + d^z 



5 . a . 



d 



dx^ 



dx^ dx^ 



(Note that these quantities differ to first order in h from 
what an observer in the lab system would measure.) One 
then obtains the following set of equations governing the 
plasma 

mn {dtv + (v ■ V)v) =jxB-Vp + g (11) 

E + vxB = (12) 



dtn + V ■ (nv) = 

Vp = kBT'Vn 
V X B = /^oj + iE 
9tB + V X E = -JB 



(13) 
(14) 
(15) 
(16) 



where 



g = -G-f gpic 



Scm — coupling 



G= T 



(w^ -c){h+Vx + hy,Vy), 



(^2 -c){hy,Vx - h + Vy), 

gprcssurc = - [h+^xP + hydyP, h^d^p - h+dyP, 0] 

gem — coupling ^ [ i^xBzJx ^ [h'+By — llxBxjJzi 

2h+Bzjx + h^Bzjy, -{h+Bx + hy,By)jy\ 
J£ = [dy{h+Bz) + dz{h+By - h^B,), 

-dx{h+Bz), -dx{h+By - h^By)] 
JB = [dt{h+B^ + hxBy),0,-dt{h+Bz)] 

We want to point out that terms of the order Eh^/c 
\n '}e have been neglected, since E/B is of the order Ca 
in ideal MHD theory. Furthermore, note that cF\ ^ 
cF\ = Ex, etc., which is the origin of some of the terms 
appearing in the expressions for gom-coupUng, 'iE and js. 
In addition to eqs.(pT|)-(pj|) Maxwell's equations produce 
constraints (e.g. for V • B), however it is easy to verify 
that these constraints are propagated by the equations 
of time evolution, (ITll). (fl and (|l6|). 



III. WAVE- WAVE INTERACTIONS 

From now on we will consider resonant three-wave 
interaction between gravitational radiation and MHD- 
waves. Such processes may occur whenever quadratic 
nonlinear terms, such as the right hand sides of eqs. ([ll|)- 
([iq), are present, and the dispersion relations of the in- 
teracting waves allow for the frequency and wave-vector 
matching conditions to be fulfilled. For a general review 
of resonant three wave interaction in plasmas, see e.g. 
Ref. 22. 

In the absence of any waves we assume to have the 
configuration of a static homogeneous, n = n'*^' , magne- 
tized, B = B^''^ plasma in Minkowski space. Cartesian 
coordinates are chosen {x^^ — {ct, x, y, z)) for a frame in 
which the velocity field (and the current density field) 
vanishes. The gravitational waves are then inferred as 
small perturbations to the Minkowski background, as in 
the previous section, and the MHD-waves as the existence 
of the smah fluctuations: n^^), v^^), j(i),E(i', B^^). Fur- 
thermore, in order to simplify the algebra, we make the 
assumptions that the direction of B'^^ is everywhere par- 
allel to the direction of propagation of the gravitational 



waves, i.e. B^") = B^^^z, and that the gravitational ra- 
diation is polarized such that /i+ = 0. 



A. Linear Calculations 

It is instructive to first investigate the linearized the- 
ory in some detail. Linearizing the equations (|ll|)-(|lq) in 
the variables h^ , n'^', v^^), j(i), E^^), B'^' we find that the 
gravitational waves do not drive plasma perturbations 
linearly. This is a consequence of the direction of prop- 
agation of the gravitational wave (parallel to the mag- 
netic field) that was chosen. Similarly the linear plasma 
perturbations are the ordinary MHD-modes. Fourier an- 
alyzing we obtain the dispersion relations for the shear 
Alfven wave 

Da^oj"^- C\kl = (17) 

and for the fast and slow magnetosonic wave 



D,n = UJ^ 



u^k^{Cl+C\) + klk^ClC\^Q. (18) 



= n(i) 



XV: 



(1) 



lWfc„ 



^z;li) 



UJr, 



^niz^AX q(1) I kmxC aX q{1) 



B(o) 



UJ„ 



B(o)l 



(22) 



for the magnetosonic modes, with the frequency-wave 
number pairs (luat'^a) and {corm'k-m) satisfying the dis- 
persion relations Da = and Dm = 0, respectively. The 
constant x is defined as x = nQ{ujm — Cgkm z) / (^'s^mx^m- 
With aid of the corresponding eigenvectors and the rela- 
tion vV'' = vi^^Clkmxkmz/i^m " ^gkl^J we Can, after 
some algebraic manipulations, write the normal modes 
as 



aA 






Ti(°) dD^ 



„(1) 



(23) 
(24) 



For the nonlinear calculation we need the eigenvectors 
expressed in terms of the normal modes, and the final 
linear results are 



The constants introduced are the Alfven velocity C\ = 

^(0) /mn^^^fio, the thermal velocity C| = ksT /m^ and 
we have used the notation k = fca;X -f fc^z together with 
k = |k|. In the next subsection we will consider superpo- 
sition of MHD- waves, and by expressing all variables in 
terms of the fluid velocity we can represent the solution 
as a sum of eigenvectors 



v(i) 


/ 


E(i) 




Vbw y 


\ 






\ 



^k„ X B(o) - ^2^k„ 



Mo 

B(o) X vi^) 



X V, 



(1) 



0„B(o) 



,(1) 



(19) 



,(1) 



/wq,, (Tc = kazB^^> lua and a is a 



where 0q = k^ • Vq 
wave-mode index. 

As we intend to study nonlinear wave coupling it is con- 
venient to adopt the normal mode method of approach 
[22], which typically simplifies the algebra in the nonlin- 
ear stage of the calculations. We define a normal mode 
as a linear combination, a^, of the dynamical quantities 
that to linear order satisfies 



and 
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Vbw / 



/ 



aA 





(0,i,0) 

7 ^1 



\ 



^ 2u:AtJ.a ' "' ZwA/Jo 

(-^,0,0) 
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where Cm 

UJm/n''°'>kn 



Cikt 



(C,o,c^ 



-Cgfc^ 
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, ,2 C'^h.2 



V(- 

/2(^^ 



re<'')fe„;^Mo 



i(0)fe 



(25) 



(26) 
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^mKnz^A^s) ^^^'^ C 



B. Nonlinear Calculations 



dtaa + ILOaO-a = 



(20) 



The dynamical quantities are now only assumed to have 
harmonic spatial dependence, i.e. V = ike,. From eq. 
(pO) the proper linear combinations are 



aA 



V, 



(1) 



LUA 



Bi}^ 



y kAzB(o) y 

for the Alfven mode, and 



(21) 



The aim of this section is to investigate the lowest order 
nonlinear influence of the gravitational radiation on the 
MHD modes described above. In particular we are inter- 
ested in the threshold value (for parametric excitation) 
of the gravitational amplitude, and the growth rates of 
the excited MHD waves. We will again assume that the 
wave vectors lies in 2:z-plane, i.e. k = fc^jX + fc^z for the 
MHD-waves, but in contrast to the case of linear wave 



modes this is a restriction made in order to simplify the 
algebra [23]. 

We consider coherent three-wave interactions, the 
three waves being one gravitational wave and two MHD- 
waves, with the matching conditions 



: W/ + UJII 



i-II 



(27) 
(28) 



where / and // arc indexing the MHD- waves. In the 
nonlinear regime the normal modes does no longer satisfy 
eq.(|2y), but rather 

dtUa + iuJaUa = {[dtaa]n.l\ 



Another special case of particular interest is the limit of 
parallel propagation (but with arbitrary ratio Cg/Cj^), 
in which case the magnetosonic dispersion relation coin- 
cides with that of the shear Alfven wave, and the only 
distinction between the modes is the polarization, which 
differ by 90 degrees. Again the general coupling coeffi- 
cient (p^) seem to diverge, since from the magnetosonic 
dispersion relation uj — C\k'j ~ fc^j when k^j — *■ 0. How- 
ever, by using another renormalization 



ai 
Ci 
Cii 



ai/x 
Ci/x 
xCn 



Ci — Cii — — 7:w„ 



where n.l. denotes (first order) nonhnear terms and the and taking the limit fc/^,, ku^ -^ 0, we obtain 
suffix ka indicates that terms not oscillating as e^^'^ van- 
ishes due to rapid oscillations. Explicit forms for the right 
hand side is found by using the original expressions for 
the normal modes - eqs. (^l]) or (|2^ ) - together with eqs. 

0-®- 

We let index / denote the magnetosonic wave pertur- 
bation, index // the Alfven wave perturbation and we 
use a complex representation (i.e. letting / ^ / + /* for 
all variables, where the star denotes complex conjugate). 
Making use of the linear eigenvectors (|2^) and (|2^) as ap- 
proximations in the nonlinear right hand sides, we obtain 
the coupled mode equations 



(35) 



dtai + iujiai = Cia*jjhy 
dto-ii + iojiiaii = Ciia*jhy 



(29) 
(30) 



after lengthy but straightforward algebra, where the cou- 
pling coefficients are [24] 



Ci^ 
Cii = 



i nokxi 
'2u;j-Clkf'' 

i 1 ujjCgkxi 

'2m{uj-ClC\k]klj) 



iVllLOg 



(31) 
(32) 



In deriving the expressions for C/ and Cn we have also 



— k//, which follows from the matching 
together with Ca ^ c and the dispersion 



applied k/ 

condition 

relations. 

From eq. (pL) and ( p2| ) one may get the incorrect im- 
pression that the coupling strength diverges in the limit 
Cg -^ 0. Thus in order to shed some light on our formu- 
las in the cold limit, we first renormalizc 



C|a/ 



Since dissipation of the waves have not been included in 
our model, the instability threshold value of the grav- 
itational amplitude found from (|2^) and ( PQ ) is so far 
zero. However, since only weak dissipation is of interest 
we can take such effects into account by simply substi- 
tuting dttta — > {dt + "fa)0'a in the coupled mode equations 
[25], where ja is the linear damping rate of the mode a. 
The most common damping mechanism of MHD waves 
is that due to finite resistivity. Calculating the linear 
damping by replacing (|lj) with E + v x B = 77J, where 
T] is the resistivity, we find 7^ — rjk'^/fio- Next we intro- 
duce the (weakly time dependent) normal mode ampli- 
tudes, Aa, defined by «„ = AqC"*""*, where a = I, II. 
Substituting these expressions into (p3) and (Bfl) tak- 
ing the damping into account, we find the general form 
for the condition of parametric growth of waves \hx\ > 
hthr = (7/7///C/C'|'j)^/^[22], where fix is the amplitude 
of the gravitational wave and hthr is the threshold value 
for parametric excitation. In the limit of parallel prop- 
agation we find from ( p5| ) that the threshold value /ithr 
reduces to 



|/ixl > h 



thr 



47 



I7II 






(36) 



ai -^ 
Ci -^ 
Cii ^ 

and then take the limit C| 
coupling coefficients then becomes 



C^sCi ^ 
Cii/Cl 



Furthermore, if the gravitational amplitude is well above 
threshold {\hx\ ^ /ithr) the general expression for the 
parametric growth rate F from (^9|) and (30) is F « 
^J(JJCYi\hy.\ [22], and the result for the special case of 
parallel propagation is 



0. The corresponding 



r« 2^g\hx\ 



(37) 
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Cii^ 
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2 noLOj 



(33) 
(34) 



It should be pointed out that in addition to the wave 
interactions considered above, we have found zero cou- 
pling coefficients for a number of cases. To be specific: 
For the same polarization of the gravitational pump wave 



(/i+ = 0), and propagation parallel to the external mag- 
netic field, the following combinations of MHD-waves 
cannot be excited in the resonant three wave approxi- 
mation, since the coupling coefhcient then becomes zero: 
1) Two ion-acoustic (or slow magnetosonic) modes. 2) 
One ion-acoustic and one Alfven wave. 3) Two Alfven 
waves with the same linear polarization. 

Note that for the case of non-zero coupling considered 
above, the Alfven waves have perpendicular polarizations 
in the parallel limit. The coupling between the differ- 
ently polarized modes then results from the quadratic 
nonlinear terms proportional to hxVyii. and h^VxY in 6q. 
(|l l|) , and from similar cross terms in (|l6| ) that couple 
the two different polarizations through hy^. The depen- 
dence of the results on the various polarizations for par- 
allel propagation can be physically understood as follows: 
For the MHD-waves to gain energy from the gravitational 
wave, the MHD-waves must be able to reduce the grav- 
itational wave amplitude. Including the source term in 
the wave equation for the gravitational wave, we immedi- 
ately see that it is only the component T^y of the energy 
momentum tensor that may affect the gravitational wave 
with /ix -polarization. However, for the three "null-cases" 
listed above, it is trivial to see that T^y vanishes (within 
the quadratic MHD approximation), and thus the gravi- 
tational wave is unaffected by the presence of such waves. 
From energy conservation [26] it is thus clear that the 
corresponding coupling must disappear. For perpendicu- 
lar polarization of the Alfven waves, on the other hand, 
Txy = mnoVxVy — B^By/ jjlq ^ 0, and thus the gravi- 
tational amplitude is affected by such a combination of 
waves and in accordance with this we have found the 
coupling to be non-zero. 



IV. SUMMARY AND DISCUSSION 

We have considered parametric excitation of Alfven 
waves by gravitational radiation propagating parallel to 
the external magnetic field. As a starting point, stan- 
dard ideal MHD equations (i.e. without special relativis- 
tic effects) incorporating the curvature of space time have 
been derived. It should be pointed out that the system 
of equations (0) -(H) in principle can be used in situa- 
tions where we have strong deviation from Minkowski 
space time, although the condition of non-relativistic 
fluid velocities then limits the applicability. Focusing on 
the case where the metric is that of a small amplitude 
monochromatic gravitational wave superimposed on flat 
space-time, the growth rate for nonlinearly coupled shear 
Alfven and fast magnetosonic waves have been found. 

In our calculations we have considered a monochro- 
matic gravitational pump wave, which could be produced 
by binary systems. As seen from (B^ (or more generally 
from (p^) and (p2) which applies for arbitrary directions 
of propagations [27]), the growth rate is roughly of the 



order F ~ h^LOg. Thus the plasma parameters no, Bq 
and T do not significantly influence the growth rate, at 
least not as long as the assumptions of the derivation is 
fulfilled. It is not hard to find a plasma fulfilling these 
assumptions, e.g. choosing no ^ 10"'^^ m~^[28] and con- 
sidering waves with frequency w < 10"* rad s^^, -Bo rnay 
attain any value roughly in the interval 10~^ — IT for 
the estimate F ~ /ix^g to apply, where the limits of the 
interval comes from the conditions uj <^ ujc and C\ ^ c^ 
respectively. Furthermore, the condition C| ^ c^ still al- 
lows for comparatively high plasma temperatures, since 
it is obviously fulfilled provided the thermal velocities of 
the particles are much smaller than the speed of light. 
As for the source of gravitational radiation we follow 
Ref. 12 and consider a binary system of two equal masses 
m — 3Mq separated by a distance of six Schwarzschild 
radii r — 12Gm/c^ emitting gravitational radiation with 
frequency of the order ujg « 10^ rad/s. By considering 
two point masses separated by a fixed distance r the am- 
plitude hx Sit a distance R from the system is estimated 
to give /ix ^ Gmr^uj^/{2c^R) . Also in the case of Ref. 12, 
which considered parametric excitation of high-frequency 
plasma waves, the growth rate fulfilled F '^ h^LOg, imply- 
ing the growth rate F ^ 10^'^ s^^ at a distance of 1/60 
au from the source, where a process at a closer distance 
was ruled out by the frequency matching conditions com- 
bined with the linear dispersion relations. In our case the 
linear dispersion relations and matching conditions allow 
a parametric process closer to the source, and thereby 
opens up the possibility for a higher growth rate, al- 
though too close to the source the background plasma 
may be too inhomogeneous and too far from steady state 
for our calculations to be applicable. 

Furthermore, excitation of MHD waves may take place 
in a dense plasma, and therefore processes such as super- 
novas are of interest, where gravitational wave absorption 
may take place inside the exploding star. In a discussion 
of possible mechanisms of absorbing gravitational wave 
energy in supernovas Ref. 14 has written "Since the ef- 
fect of acceleration by gravitational waves is independent 
of mass of the charge, both the ions and the electron re- 
spond in an identical manner, which is not the case for 
electromagnetic waves. This means that waves such as 
Alfven waves which describe oscillation of charge neu- 
tral plasmas are ideal. The coupling, however, is weak." 
At the present stage of understanding it is too early to 
deduce whether significant gravitational wave absorption 
by MHD waves may occur. Calculations taking into ac- 
count the effects of a broad band gravitational spectrum, 
plasma inhomogeneities, etc., must first be performed. 
In particular inhomogeneity scale lengths with a scale 
length significantly shorter than the wavelength of the 
gravitational mode - such as at the plasma boundary 
of the supernova - may lead to excitation of MHD sur- 
face waves with a significantly enhanced growth rate as 
compared to the present homogeneous plasma coupling 



mechanism. This is in analogy with parametric excita- 
tion scenarios for high frequency plasma surface waves 
[29], where the surface waves may have a considerably 
higher growth rate than the corresponding bulk waves, 
provided the inhomogeneity scale length is considerably 
shorter than the wave length of the pump wave. Such a 
problem, however, is a project for future work. 

An interesting question from a theoretical point of 
view is whether the coupling coefficients Cj and C// of 
eqs. (Bl]) and (p2h satisfy the Manley-Rowe relations 
[22]. Generally such relations follow from an underlaying 
Hamiltonian structure of the governing equations, and 
assures that each of the decay products takes energy from 
the pump wave in direct proportion to their respective 
frequencies. This means that the parametric process can 
be interpreted quantum mechanically - i.e. we can think 
of a three wave process as the decay of a pump wave 
quanta with energy hujg into wave quantas with energy 
Tiwi and hujji respectively. An interesting consequence 
is that generally a lot of three wave decay processes are 
forbidden from the start by the Manley-Rowe relations 
(for example the decay of a plasmon into two photons), 
since they imply that we only get a positive growth rate 
when the pump wave has the highest frequency, in consis- 
tence with the quantum picture [30] . It should be pointed 
out that all well established basic systems of equations 
in plasma physics such as the ideal MHD-equations, the 
Vlasov-Maxwell equations and the standard multi-fluid 
equations, all possess the underlaying Hamiltonian struc- 
ture that leads to fulfillment of the Manley-Rowe rela- 
tions for three- wave interaction processes [31]. However, 
since the theory of gravitation differs in important re- 
spects from other fundamental theories of physics, it is 
an open question whether the Manley-Rowe relations ap- 
plies also when one of the interacting modes is a gravi- 
tational wave. To investigate this issue we consider the 
energy increase rate of mode I. Multiplying Eq. (|29| ) with 
-k1jC\)/L0ikxi and using the eigenvectors 



(|2q) we find 

dWi 

dt 



ujiV 



(38) 



where V = 2lTa[mnQijjghy^v*^^v*jj/uJi\ and the mag- 

mnQk'j(ujj — 

cf. Eq. (20) in Ref 

with mnQV*jj, and 



netosonic wave energy density is Wj 
kljC\){dD^/duJi) \v.,i\^ I2kl,uji 
32). Similarly, multiplying Eq. ( 
applying (E5h and (EQ) we obtain 



dW, 



II 



dt 



= ujiiV 



(39) 



where Wn = mriQ \vyii\ (cf. Eq. (13) in Ref 32). Up to 
now, we have treated the gravitational wave as a pump 
wave, i.e. we have only considered the initial stage of 
an instability where the energy density of the daughter 
waves are small enough for the influence on the gravita- 
tional wave to be neglected. For a practical purpose this 



regime may apply until nonlinear saturation mechanisms 
(outside our calculation scheme) sets in for the growing 
MHD-waves, unless the plasma is extremely dense. From 
a theoretical point of view, however, it is of interest to 
consider the decrease in the gravitational wave amplitude 
due to the growth of the MHD perturbations. Including 
the energy momentum source term in the wave equation 
for the gravitational wave, i.e. using 'Oih^i, = —2kT^i, = 
— IGttGT^j^/c'*, keeping only the resonantly varying part 



of T^i, (i.e. the part proportional to exp[i(fct,z 



M) 



and letting hi2 = /ix (i) exp[i(fcgZ — ojgt)] + c.c, we im- 
mediately obtain 



IWo 



9/ix 
dt 



k{T, 



12jk„ 



mnaVxiVyii 



BxiByll 



A^o 



(40) 



which, after using the eigenvectors (Eq) and (Eq) and mul- 
tiplying with Wg/i* can be written 



dt 



-U^gV 



(41) 



where Wg — uj'^\hx 



"^ /2k is the energy density of the grav- 
itational wave [33]. Since the same factor V appears in 
Eqs (H), (H) and (|l|), the Manley-Rowe relations are 
indeed fulflUcd (i.e. each mode changes energy in direct 
proportion to its frequency) , and furthermore we see that 
the total wave energy W = Wj + Wu + Wg is conserved, 
which follows from (|3g), ( |39| ) and (|4|) together with the 
frequency matching condition ([2 7]). 
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